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Abstract 

We formulate a statistical model for description of nuclear composition and 
equation of state of stellar matter at subnuclear densities and temperature up to 20 
MeV, which are expected during the collapse and explosion of massive stars. The 
model includes nuclear, electromagnetic and weak interactions between all kinds of 
particles, under condition of statistical equilibrium. We emphasize importance of 
realistic description of the nuclear composition for understanding stellar dynamics 
and nucleosynthesis. It is demonstrated that the experience accumulated in studies 
of nuclear multifragmentation reactions can be used for better modelling properties 
of stellar medium. 

PACS: 26.50.+X , 21.65.-f, 25.70.Pq , 26.30.-k, 97.60.Bw 
1. Introduction 

In violent nuclear reactions strong interaction between many nucleons leads to a fast 
equilibration. This short-range interaction is responsible for a sharp freeze-out when 
the inter-particle distance becomes larger than the interaction range. For these reasons 
statistical models have proved to be very successful for interpretation of nuclear reactions 
at various energies. They are widely used for description of fragment production when one 
or several equilibrated sources can be identified. Originally this concept was proposed for 
compound nucleus decays, such as evaporation or fission of excited nuclei IJJ. Recently, 
it was demonstrated that the concept of equilibrated source can even be effectively used 
for more violent multifragmentation reactions leading to production of many fragments 
p] . On other side, the statistical equilibrium is expected in many astrophysical processes, 
when the characteristic time for formation of nuclei is much shorter than the time-scale 
of these processes. For example, one of the most spectacular astrophysical events is a 
type II supernova explosion, with a huge energy release of about several tens of MeV per 
nucleon [Sill]. When the core of a massive star collapses, it reaches densities several times 
larger than the normal nuclear density po = 0.15 fm~'^. The repulsive nucleon-nucleon 
interaction gives rise to a bounce-off and creation of a shock wave propagating through the 
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in-falling stellar material. This shock wave is responsible for the ejection of a star envelope 
that is observed as a supernova explosion. During the collapse and subsequent explosion 
the temperatures T ^ (0.5 -7- 10) MeV and densities p ~ (10~^ -r- 2)po can be reached. It 
is widely believed that the nuclear statistical equilibrium should be reached under these 
conditions. As shown by many theoretical studies, a liquid-gas phase transition should 
take place in nuclear matter under such conditions. 

As discussed by several authors (see e.g. [SI El El E] ) , there are problems in producing 
successful explosions in hydro dynamical simulations of the core-collapse supernovae, even 
when neutrino heating and convection effects are included. The hope is that full 3-0? 
simulations will help to solve this problem [9]. On the other hand, it is known that 
nuclear composition is extremely important for understanding the physics of supernova 
explosions. In particular, the weak reaction rates and energy spectra of emitted neutrinos 
are very sensitive to the presence of heavy nuclei (see e.g. [IHl [TH [T21 US])- This is also 
true for the equation of state (EOS) used in hydrodynamical simulations since the shock 
strength is diminished by dissociation of heavy nuclei. 

Nuclear reactions in a supernova environment are especially important because super- 
nova explosions may be considered as breeders for creating chemical elements. According 
to present understanding, there are three main sources of chemical element production 
in the Universe. The lightest elements (up to He, and, partly, Li) are formed during the 
first moments of the Universe expansion, immediately after the Big Bang. Light elements 
up to ^^O can be produced in thermonuclear reactions in ordinary stars like our Sun, 
while heavier elements up to Fe and Ni can be formed in heavy stars at the end of the 
nuclear burning epoch. It is most likely that heavy elements up to U were synthesized in 
the course of supernova explosions. Pronounced peaks in the element abundances can be 
explained by neutron capture reactions in s- and r-processes IT5]. It is believed that 
suitable conditions for the r-process were provided by free neutrons abundantly produced 
in supernova environments together with appropriate seed nuclei. 

The EOS of supernova matter is under investigation for more than 25 years. One of the 
first EOS, frequently used in supernova simulations, was obtained in refs. (THl [IT] many 
years ago. It includes both light and heavy nuclei in statistical equilibrium. However, 
it does not include the whole ensemble of hot heavy nuclei, replacing them by a single 
"average" nucleus. The same assumption within a relativistic mean-field approach was 
used in the EOS of ref. [IH]. As was already pointed out by many authors [IHl [201 1211 122] 
this assumption is not sufficient for accurate treatment of the supernova processes. We 
think that this kind of approximation can distort the true statistical ensemble in many 
cases. There are other statistical calculations which consider the ensemble with different 
nuclear species, but in the partition sum they include only nuclei in long-lived states 
known from terrestrial experiments (see refs. [ini[2Il[22])- Also, for description of unknown 
neutron- rich hot nuclei only properties (e.g., the symmetry energies) of cold and slightly 
excited isolated nuclei have been used up to now. These assumptions are not justified for 
supernova environments characterized by relatively high temperatures (up to 10 MeV) 
and densities of electrons and baryons in the range p ^ 10~^ -7- 10~^po- We believe, 
in order to achieve a more realistic description of supernova matter, it is necessary to 
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use rich experience accumulated in recent years by the nuclear community in studying 
highly excited equilibrated systems in nuclear reactions. In particular, multifragmentation 
reactions provide valuable information about hot nuclei in dense surrounding of nucleons 
and other nuclei, which in many aspects is similar to supernova interior [23] . 

2. Nuclear reactions in supernova environments. 

2.1 Studying equilibrated nuclear matter in laboratory. 

Properties of strongly-interacting nuclear matter are studied experimentally and theo- 
retically for a long time. At present there exist consensus about phase diagram of nuclear 
and neutron matter (see, e.g., refs. |2H [25l |T6|). It is shown schematically in Fig. 1 
for symmetric nuclear matter, for range of densities and temperatures expected in the 
Supernova II explosions. It is commonly accepted that this diagram contains a liquid- 
gas phase transition. From this phase diagram one can conclude that nuclear matter 
at densities p ~ 0.3 — 0.8po and temperatures T < Tc should be in the mixed phase. 
This phase is strongly inhomogeneous with intermittent dense and dilute regions. In the 
case of electrically-neutral matter this mixed phase may have different topologies such 
as spherical droplets, cylindrical nuclei, slab-like configurations and others. These con- 
figurations are generally referred to as nuclear 'pasta' phases [2B|, which were recently 
under intensive theoretical investigation [271 [28l [29]. However, in the coexistence region 
at lower densities, p < O.Spo, which are considered in this paper, the nuclear matter 
breaks up into compact nuclear droplets surrounded by nucleons. These relatively low 
densities dominate during the main stages of stellar collapse and explosion. Description 
of nuclear composition in this region requires theoretical extrapolation of nuclear proper- 
ties to these extreme conditions. As became obvious after intensive experimental studies 
of nuclear multifragmentation reactions, they proceed through formation of thermalized 
nuclear systems characterized by subnuclear densities p ~ O.lpo and temperatures of 3-8 
MeV. Thermodynamic conditions associated with these reactions are indicated by the 
shaded area in Fig. 1. This gives us a chance to extract properties of hot nuclei in the 
environment of other nuclear species directly from the experimental multifragmentation 
data and then use this information for more realistic calculations of nuclear composition 
of stellar matter. We have also shown isentropic trajectories with the entropy per baryon 
(S/B) of 1, 2, and 4 units (see section 5.2) typical for supernova explosions. One can 
see, for example, that an adiabatic collapse (and expansion) of stellar matter with typical 
entropies of 1-2 units per baryon passes exactly through the multifragmentation domain. 

Nuclear multifragmentation, i.e. break-up of hot heavy nuclei into many fragments, 
was under intensive investigation during the last 20 years. It was solidly established by 
both theoretical and experimental studies that this channel dominates at high excitation 
energies, above 3-4 MeV per nucleon, replacing sequential evaporation and fission of the 
compound nucleus, which is conventional mechanism at low excitation energies. In this 
respect, multifragmentation is a universal process expected in all types of nuclear reac- 
tions, induced by hadrons, heavy ions, and electromagnetic collisions, where the nucleus 



3 



receives a high excitation energy. The Statistical Multifragmentation Model (SMM) [2] 
is one of the most successful models used for theoretical interpretation of these reactions. 
Some examples, how low-density equilibrated nuclear systems can be produced, and how 
well their decay can be described within this statistical approach, can be found in refs. 
[2l[3Ql[3ll|32l[Sll3ll35l[36]. Recently, experimental evidences have appeared [371 EHl EH] 
that the symmetry energy of hot fragments produced in multifragmentation reactions is 
significantly reduced as compared with the value in cold nuclei. The surface and bulk 
energy of nuclei in this hot and dense environment can be modified too [101 SI] • These 
conclusions, made after appropriate analyses of experimental data, suggest that the same 
modifications will occur in stellar matter at similar densities and temperatures. Actually, 
the "in- medium" modifications are quite expected, since the nuclei will interact with the 
surrounding matter, and, therefore, should change their properties. These modifications 
may have important consequences for nuclear composition, equation of state, and weak 
reaction rates on these nuclei. 

2.2 Nuclear and electro-weak reaction rates. 

In the supernova environment, as compared to the nuclear reactions, several new 
important ingredients should be taken into consideration. First, the matter at stellar 
scales must be electrically neutral and therefore electrons should be included to balance 
a positive nuclear charge. Second, energetic photons present in hot matter may change 
nuclear composition via photo-nuclear reactions. And third, the matter is irradiated by 
a strong neutrino wind from the protoneutron star. 

We consider macroscopic volumes of matter consisting of various nuclear species with 
mass number A and charge Z, {A, Z), nucleons {n = (1, 0) and p = (1, 1)), electrons (e~) 
and positrons (e^) under condition of electric neutrality. There exist several reaction 
types responsible for the chemical composition in supernova matter. At low densities and 
temperatures around a few MeV the most important ones are: 1) neutron capture and 
photodisintegration of nuclei, which proceed via production of a hot compound nucleus 

{A, Z)+n^{A + l,Zy ^{A + l,Z)+j , ... 

{A,Z)+^^iA,Zr^iA~l,Z)+n, ... (1) 

2) neutron and light charged particle emission (evaporation) by the hot nuclei 

{A,Zy^{A-l,Z)+n, {A,Zy^{A-l,Z-l)+p, ... (2) 

and 3) weak processes induced by electrons/positrons and neutrinos/antineutrinos 

{A, Z)+e- ^ {A,Z + u , {A, Z) + e+ ^ {A, Z + 1) + i> , (3) 

which transfer protons to neutrons and vice versa. There are many other reactions not 
shown here, which are naturally taken into account within the assumption of statistical 
equilibrium. The characteristic reaction times for neutron capture, photodisintegration 
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of nuclei and nucleon emission are defined as 



Tcap = [{(^nAVnA) Pn]~ 
T<yA = [{(^■yAVjA)P'y] 



-1 



-1 



n,p ; 



(4) 



respectively. Here anA and cr^^ are the corresponding cross sections, VnA and v^a are the 
relative (invariant) velocities, and Tn,p is the neutron (proton) decay width. 

In our calculations for anA we use the geometrical neutron-nuclear cross sections, that 
is a good approximation for the considered range of temperatures (T ^ (0.5 ^ 10) MeV). 
The photo-nucleus cross section a^A was taken phenomenologically under assumption 
that it is dominated by the giant dipole resonance. These parametrizations of neutron 
and photon cross-sections are in a good agreement with experimental data (see discussion 
in ref. [2]). The evaporation decay widths were calculated according to the Weisskopf 
evaporation model as described in ref. [12] . Our estimates show that at temperatures and 
densities of interest these reaction times vary within the range from 10 to 10^ fm/c, that 
is indeed very short time scale compared to the characteristic hydrodynamic time of a 
supernova explosion, about 100 ms [3]. The nuclear statistical equilibrium is a reasonable 
approximation under these conditions. 

We have calculated the reaction rates of Eq. (jl]) for nuclei with A =60 and Z =24, 
which are typical for stellar nucleosynthesis in dense matter and at a typical electron 
fraction (i.e., the ratio of electron and baryon densities Pe/ps) Ye ~ 0.4. They are 
presented in Fig. 2 as function of neutron density for several temperatures. By analyzing 
this figure one should take into account that the neutron density p„ is usually by 2-5 
times smaller than p. One can see clearly that for densities p„ > 10~^po and for the 
expected temperatures of the environment, T<5 MeV, we obtain t^a » Tcap,Tn,p, i-e. 
the photodisintegration is more slow than other processes. There exists a range of densities 
and temperatures, for example, pn > 10~^po at T = 1 MeV, pn ^ 10~^po at T = 3 MeV, 
and p„>10~^po at T = 5 MeV, where the neutron capture dominates, i.e. Xcap < Tn^p. 
Under these conditions new channels for production and decay of nuclei will appear (e.g. 
a fast break-up with emission of a-particles or heavier clusters) which restore the detailed 
balance. We expect that in this situation an ensemble of various nuclear species will 
be in chemical equilibrium like in a liquid-gas coexistence region, as also observed in 
the multifragmentation reactions. Here the nuclear system is fully characterized by the 
temperature T, density p (which is nearly the same as baryon density pb), and electron 
fraction Ye. One may expect that modifications of nuclear properties come into force in this 
environment, because of intensive interaction between clusters. This is complementary 
to the well known effects in isolated nuclei: at high temperature the masses and level 
structure in hot nuclei can be different from those observed in cold nuclei (see, e.g.. 



The weak interaction reactions are much slower. The direct and inverse reactions in 
Eq. ([3]) involve both free nucleons and all nuclei present in the matter. It is most likely 
that at early stages of a supernova explosion neutrinos/antineutrinos are trapped inside 
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the neutrinosphere around a protoneutron star [H]. In this case we should impose the 
lepton number conservation condition by fixing the lepton fraction Y^. Then one should 
take into account the continuous neutrino fiux out of the surface of the neutrinosphere 
propagating through the hot bubble. Due to large uncertainties in the weak interaction 
rates, below we consider three physically distinctive situations: 

1) fixed lepton fraction Yl corresponding to a /3-equilibrium with trapped neutrinos 
inside the neutrinosphere (early stage of the explosion); 

2) fixed electron fraction Y^ but no /3-equilibrium inside a hot bubble (early and 
intermediate times); 

3) full /^-equilibrium without neutrino (the late times of the explosion, after neutrino 
escape). 

The second case corresponds to a non-equilibrium situation which may take place in 
the bubble, before the electron capture becomes efficient. Actually, this case is considered 
as basic for calculations of nuclear composition in the hot bubble behind the shock. 
Generally, one should keep in mind that weak reactions are often out of equilibrium. Our 
estimates show that their characteristic times range from 10 ms to 10 s depending on 
thermodynamical conditions and intensity of the neutrino wind. Therefore, one should 
specify what kind of statistical equilibrium is expected with respect to weak interaction. 

3. Formulation of the statistical model. 

Below we describe supernova matter as a mixture of nuclear species, electrons, photons, 
and perhaps neutrinos in thermal equilibrium. For the macroscopic scales one can safely 
apply the grand-canonical approximation. We call this model the Statistical Model for 
Supernova Matter (SMSM). It was first proposed in ref . j20j . 

3.1 Equilibrium conditions. 

Within the SMSM each particle i with baryon number Bi, charge Qi and lepton number 
Li is characterized by a chemical potential yUj, which can be represented as 

/ij = BijJB + QifJ'Q + LifiL. (5) 

Here /i^, /ig and fii are three independent chemical potentials which are determined from 
the conservation of total baryon number B = J2iBi electric charge Q = Yl,iQi and lepton 
number L = YliLi of the system. Explicitly, the chemical potentials for nuclear species 
(A, Z), electrons (e~, e"*"), and neutrinos (z/, z>) can be expressed as 

/iAZ = A^J'B + Z^lQ , 

yWe- = -/ie+ = -/^Q + , (6) 
i^u = -fJ-u = IJ'L ■ 

These relations are valid also for nucleons, fin = f^B and fip = l^s + I^q- 
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The corresponding conservation laws can be written as 
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(7) 



PL = Pe + Pu - Pu = YlPB ■ 



Here pe = Pe- — Pe+ is the net electron density, is the lepton fraction. The second 
equation requires that any macroscopic volume of the star is electrically neutral. 

The lepton number conservation is a valid concept only if u and v are trapped in 
the system within the neutrinosphere [44j. If they escape freely from the system, the 
lepton number conservation is irrelevant and pi = 0. In this case two remaining chemical 
potentials are determined from the conditions of baryon number conservation and electro- 
neutrality. Since the (3— equilibrium may not be achieved in a fast explosive process, we 
also often fix the electron fraction Yg = Pe/ps in this case. 



3.2 Ensemble of nuclear species. 



Our treatment of nuclear reactions is based on the statistical multifragmentation model 
(SMM) [2], which was very successfully applied for description of experimental data. For 
describing an ensemble of nuclear species under supernova conditions we use the Grand 
Canonical version of the SMM [15] . After integrating out translational degrees of freedom 
one can represent pressure of nuclear species as 



TY,Paz = TY.9Az—^ir^^P 
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(8) 



where pAz is the density of nuclear species with mass A and charge Z. Here qaz is the 

1 /2 

ground-state degeneracy factor of species {A,Z), At = (27rhymNT) is the nucleon 
thermal wavelength, rriN ^ 939 MeV is the average nucleon mass. V is the actual volume 
of the system, and Vf is so called free volume, which accounts for the finite size of nuclear 
species. We assume that all nuclei have normal nuclear density po, so that the proper 
volume of a nucleus with mass A is A/po. At low densities the finite-size effect may be 
included via the excluded volume approximation Vj/V ~ (1 — Pb/po)- We emphasize 
that this last approximation is commonly accepted in statistical models. As we know 
from nuclear multifragmentation studies information about free volume can be extracted 
from analysis of experimental data [32]. At densities pb > O.lpo the extracted Vf may 
slightly deviate from above approximation, however, they are in qualitative agreement. 
In the present work we remain in the framework of the conventional statistical approach, 
although allowing for modifications of nuclear properties. 

The internal excitations of nuclei play an important role in regulating their abundance, 
since they increase significantly their entropy. Some authors (see, e.g., ref. [12]) limit the 
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excitation spectrum by particle-stable levels known for low excited nuclei. Within the 
SMM we follow quite different philosophy. Namely, we calculate internal excitation of 
nuclei by assuming that they have the same internal temperature as the surrounding 
medium. In this case not only particle-stable states but also particle-unstable states will 
contribute to the excitation energy and entropy. This can be justified by the dynamical 
equilibrium of nuclei in hot environment, and supported by numerous comparisons with 
experiment (see part 2.1). Moreover, in the supernova environment both the excited 
states and the binding energies of nuclei will be strongly affected by the surrounding 
matter. By this reason, we find it more appropriate to use an approach which can easily 
be generalized to include in-medium modifications. Namely, the internal free energy of 
species {A, Z) with A > 4 is parameterised in the spirit of the liquid drop model, which 
has been proved to be very successful in nuclear physics [Ij: 

F^ziT, p) = Ff^ + Fl, + FJ- + F^^ . (9) 

Here the right hand side contains, respectively, the bulk, the surface, the symmetry and 
the Coulomb terms. The first three terms are taken in the standard form [2], 

(10) 




nziT)=(3o\,f^,^ A'" , (11) 

Flz=l^-^^^P^ , (12) 

where wq = 16 MeV, Eq = 16 MeV, (3o = 18 MeV, % = 18 MeV and 7 = 25 MeV are the 
model parameters which are extracted from nuclear phenomenology and provide a good 
description of multifragmentation data [2[ [31], [321 [33l [SH |35l [36] . However, these param- 
eters, especially the symmetry coefficient 7, can be different in hot nuclei at multifrag- 
mentation conditions, and they should be determined from corresponding experimental 
observables (see discussion in refs. jlQ], [H] HH]). 

In the electrically-neutral environment the nuclear Coulomb term should be modified 
to include the screening effect of electrons. Within the Wigner-Seitz approximation with 
constant electron density it can be expressed as 

F^zip) = Icip)^, , (13) 



c(p) 




1/3 




where tq = 1.17 fm and pop = {Z/A)pq is the proton density inside the nuclei. The 
screening function c(p) is 1 at Pe = and at pe = pop. Here one can also use an 
approximation Pe/pop = Pb/Po, as in ref. [TB]. In our calculations, we have checked that 
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these two choices lead to very similar results, especially at small densities. We want to 
stress that both the reduction of the surface energy due to the finite temperature and the 
reduction of the Coulomb energy due to the finite electron density favour the formation of 
heavy nuclei. Nucleons and light nuclei {A < 4) are considered as structure-less particles 
characterized only by exact masses and proper volumes [2] . Their Coulomb interaction is 
taken into account within the same Wigner-Seitz approximation. 

As follows from Eq. ([H]), the fate of heavy nuclei depends strongly on the relationship 
between F^z and fiAz- In order to avoid an exponentially divergent contribution to the 
baryon density, at least in the thermodynamic limit {A oo), inequality F^z ~ /^az must 
hold. The equality sign here corresponds to the situation when a big (infinite) nuclear 
fragment coexists with the gas of smaller clusters [17] . When Faz > fJ'Az, only small 
clusters with nearly exponential mass spectrum are present. However, there exist a region 
of thermodynamic quantities corresponding to F^z ~ A^az when the mass distribution of 
nuclear species is a power-law A~^ with r ^ 2. The advantage of our approach is that 
we consider all the fragments present in this transition region, contrary to the previous 
calculations fiGl [T7] , which consider only one "average" nucleus characterizing the liquid 
phase. 



3.3 Electromagnetic and weak processes. 



At T, fi > rrie the pressure of the relativistic electron-positron gas can be written as 



P. 



247r2 



+ 



15 I /ie , 




(14) 



where first-order correction (~ rrig) due to the finite electron mass is included, (76=2 is the 
spin degeneracy factor. The net number density pe and entropy density can be obtained 
from standard thermodynamic relations, pe = dPe/dfie, and Se = dPe/dT, which give 
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(15) 
(16) 



Electron neutrinos and antineutrinos are taken into account in the same way, but as 
massless fermions, and with the degeneracy factor twice smaller than for the electrons, 
i.e., 5'i/=l. The photons are always close to the thermal equilibrium, and they are treated 
as massless Bose gas with zero chemical potential. The corresponding density p^, energy 
density e^, pressure P^, and entropy density of photons gas are given by standard 
formulae: 



P7 



4e^ 

3r' 



(17) 



where 9-^=2. 
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All kinds of particles (nuclei, baryons, electrons, neutrinos, photons) contribute to the 
free energy, pressure and other thermodynamical characteristics of the system, and we 
sum up all these contributions. Within the model we calculate densities of all particles 
self-consistently by taking into account the relations between their chemical potentials. 

3.4 Comparison with the Lattimer-Swesty model. 

We have performed calculations for sets of physical conditions expected during the 
collapse of massive stars and subsequent supernova explosions. We take baryon number 
B =1000 and perform calculations for all fragments with 1< A <1000 and 0< Z < A 
in a box of fixed volume V. This volume is determined by the average baryon density 
Pb = B/V . This restriction on the size of nuclear fragments is fully justified in our case, 
since fragments with larger masses {A >1000) can be produced only at very high densities 
p>0.5po IHlIS], which are appropriate for the regions deep inside the protoneutron star, 
and which are not considered here. 

In the beginning we compare our model with calculations within other models on the 
marked. It is necessary to mention that all models treat electrons and photons in the same 
way, therefore, differences appear entirely due to different description of nuclear species. 
One can expect that in the case of domination of radiation processes the results will be 
similar. 

Most supernova simulations are performed with the equation of state of Lattimer- 
Swesty (LS) [inillT], where an ensemble of heavy nuclei is replaced by a single "average" 
nucleus. Therefore, one can compare only integral characteristics of the stellar matter. In 
Fig. 3 we compare our SMSM results with the LS model. We show entropies, pressure, 
mass fractions of alpha particles {Xaipha)^ and heavy fragments {Xheavy)^ at different 
temperatures, and at a fixed electron fraction, versus densities. In the SMSM Xheavy 
includes all fragments with A >4, whereas in the LS model it is only the share of the 
"average" heavy nucleus. The LS calculations were taken from ref. [22]. One can see that 
the average thermodynamical characteristics (pressure, entropy) are very close in the two 
models. This remains also true if we extend comparison to other models, e.g., reported 
in [221 [19]. However, as seen from Fig. 3, mass fractions of nuclei are very different. As 
was mentioned in [22], a small yield of alpha clusters in the LS model may be caused by 
mistakes in calculations of the Coulomb corrections to their binding energies. We stress 
again that the details of the nuclear composition are very important for dynamics of the 
explosion, since it influences the total energy balance, and determines the weak reaction 
rates. 

4. Composition of matter. 

In this section we present the SMSM results concerning nuclear and lepton compo- 
sition of stellar matter, which are important for determining general thermodynamical 
characteristics of the matter, such as energy density, pressure, and entropy. The energy 
deposition into the matter with photons and neutrinos produced by external sources is 
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also considerably influenced by this composition. We pay special attention to the con- 
tribution of nuclear species whose properties may be modifled in dense environment as 
follows from recent findings in nuclear multifragmentation reactions. 

4.1 Electron and neutron fractions of stellar matter. 

Within the SMSM we can calculate the electron fraction in the electrically-neutral mat- 
ter under assumption of full /3-equilibrium. The appropriate astrophysical sites, where 
this may happen, are the relatively slow collapse stage, and the very late stages of the 
explosion, after cooling down the matter and neutrino escape. In this case we calculate 
self-consistently densities of all species by using relations between their chemical poten- 
tials, Eq. IQ, with fiL=0- The net electron density, which is equal to the proton density, 
is explicitly given as a function of the chemical potential in Eq. ffT^ . Fig. 4 presents the 
fractions of free electrons and neutrons as functions of baryon density. 

On the left panels the results are shown for the /3-equilibrium neutrino-less matter at 
T =1 and 3 MeV. One can clearly see two general trends: with increasing baryon density 
the electron fraction Ye gradually decreases and the neutron fraction Yn increases. At 
small densities, which correspond to low /ig, the electrically- neutral matter tends to be 
isospin symmetric, with a large amount of electrons. In the case of low temperatures 
(T < 1 MeV) the protons are captured in most bound nuclei with A ~ 50-60. As was 
realized long time ago [1], at large densities the electrons are absorbed by protons in the 
inverse /3-decay process, that is driven by high electron chemical potential. When we 
increase temperature (T>3 MeV) the nuclei dissociate into protons and neutrons, that 
helps to capture electrons at large densities also. At the same time, the number of free 
neutrons increases rapidly at higher baryon densities for both low and high temperatures. 
This is important for maintaining a high rate of nuclear reactions to generate equilibrium 
ensemble of nuclei. At low temperatures (T ~ 1 MeV) a noticeable change in the trend is 
seen below pb ~ 10~^po- At these densities most neutrons are bound in large nuclei, which 
are still present in the matter. For example, Y^ ~ 0.2 means that 80% of neutrons are 
trapped in the nuclei. At even lower densities, when heavy nuclei disappear (see Fig. 3) 
the number of free neutrons increases. The same trend at higher densities is explained by 
the fact that more and more neutrons are dripping out of nuclei, since the matter contains 
less and less protons. For example, at pb > 10"'^ po more than half of the neutrons are 
free. This behavior correlates with decreasing the number of electrons and a relatively 
small share of heavy nuclei in the system. At higher densities, the structure of matter 
may change because of the neutrino/antineutrino and electron/positron capture reactions 

The right panel of Fig. 4 shows the results obtained under the condition of lepton 
number conservation, i.e., at fixed values of lepton fractions Yl=0.1, 0.2, 0.3. These 
values are consistent with uncertainties concerning the neutrinosphere radius discussed 
in literature. One can see that in this case the number of free neutrons always drops 
with density refiecting formation of very big nuclei and transition to the liquid phase at 
Pb Po- While the electron fraction stays nearly constant, exhausting around 80-90% of 
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the total Yl. These results show that weak reactions affect significantly the composition 
of supernova matter. 

4.2 Mass fractions of light and heavy nuclei. 

As well known, at low densities and temperatures the nuclear matter exists in the form 
of isolated nuclei and nucleons. At terrestrial conditions the nuclei capture electrons and 
become atoms. However, at supernova conditions the atoms are fully ionized, therefore, 
the nuclei are embedded in more or less uniform background of electrons and neutrons. 
This surrounding to a large extent determines the nuclear composition of stellar matter. 

Figures 5 and 6 demonstrate the mass fractions of nuclear matter contained in heavy 
nuclei (with mass numbers A > 4), a-particles, neutrons and protons, for different elec- 
tron fractions Y^.. One can see that at low temperatures (T < 1 MeV) the matter is 
mainly composed of heavy nuclei. If the share of electrons is small, the free neutrons 
are also present. With increasing temperatures the heavy nuclei gradually disintegrate 
into a's, neutrons and protons. At low densities this disintegration happens already at 
moderate temperatures T ~ 1-2 MeV, while at subnuclear densities (p ~ O.lpo) the heavy 
nuclei survive even at higher temperatures, though they become very excited. One should 
bear in mind, however, that in this case we are dealing with the dynamical equilibrium 
between decay of excited nuclei and absorption of surrounding nucleons, as regulated by 
the reaction rates presented in Fig. 2. An interesting observation is that heavy nuclei first 
break-up into light clusters (like a) and then these clusters dissolve into nucleons with 
temperature. This is clear seen in the yields of a-particles, which demonstrate a 'rise and 
fair behaviour both with increasing temperature and decreasing density. As seen from 
Figs. 7 and 8, small clusters with A=2 and 3 are also produced in the transition region. 

4.3 Nuclear mass distributions. 

The properties of heavy nuclei are very important for understanding processes taking 
place in stellar matter. In the bottom panel of Fig. 7 we show mass distributions of nuclear 
species at pb = 10~^po and several temperatures. At low temperatures the distribution 
of heavy nuclei looks like a Gaussian with a well defined peak (see also ref . [20j ) . In this 
case the average thermodynamic characteristics of the system may not be much different 
from the ones calculated under assumption of an "average" nucleus as in ref. [17], see 
Fig. 3. However, even in this case, the width of the distributions may be important 
for calculations of weak reactions in matter. By increasing temperature we move into 
a coexistence region of the nuclear liquid-gas phase transition: The mass distributions 
become 'U-shaped', and they contain all nuclei from light to heavy ones. At higher 
temperatures the mass distributions have exponential shape. These distributions can not 
be even approximately characterized by an "average" nucleus. This evolution of fragment 
mass distributions is well established in nuclear multifragmentation reactions [21 EE]- 

The average charge-to-mass ratio for all nuclei is demonstrated in the top panel of 
Fig. 7. For this quantity both the symmetry and Coulomb energies of fragments are 
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crucially important. The charges of heavy fragments show rather regular behaviour: The 
Zj A ratio decreases slowly with mass number because of Coulomb interaction. It decreases 
less rapidly than in multifragmentation reactions due to the screening effect of electrons. 
From our calculations we came to the conclusion that the charge distribution of nuclei at 
given A can be approximated by a Gaussian, with the width az ~ \J AT /9>'~^ [l2l HSl HE] , 
where 7 is the coefficient in the symmetry energy (see Eq. fll2p ). 

We have found that the phase transition from heavy nuclei to light fragments (nuclear 
gas) always proceeds through the same sequence of mass distributions: 'U-shape', power- 
law, and exponential ones, both with increasing temperature and decreasing density. This 
opens the possibility to study the critical behaviour in stellar matter, in the same way 
as was previously done in multifragmentation reactions [3^ One can find examples 
of the temperature-driven transitions in refs. [201 [23] and in Fig. 7 (see also Fig. 22). In 
Fig. 8 we demonstrate an example of the density- driven transition at fixed temperature 
of 1 MeV. As we can clearly see, the mass distributions evolve from 'U-shape' at densities 
(in units lO^^po) 1-0 and 0.32 to power-law at 0.18 and 0.1, and finally to exponential at 
0.03. At density 0.18 the distribution of large clusters is most fiat, that may be considered 
critical point of the phase transition [TF] . 



4.4 Shell effects. 



Up to now we have used the liquid-drop description of nuclei in the stellar matter. 
However, it is well known that at low temperatures (T < 1 MeV) the shell corrections to 
nuclear masses becomes important, at least, at terrestrial densities of matter. It is likely 
that the shell effects may also play a role in stellar matter at subnuclear densities. In this 
section we demonstrate how shell effects may influence the mass and charge distributions 
of nuclei. In particular, we have analyzed possible existence of superheavy elements in 
stellar matter, assuming that they are sufficiently long-lived. It was assumed that there 
is a certain shell correction to the free energy Faz of a specific nucleus, AF^^, taken as 
follows: 



AZ 



(iV-iV,,)2\ ^ ^ ( {Z-Z,r. 



Eshcxp + Eshcxp 



;i8) 



We have assumed existence of the island of stability around Zs/j=120 and A^sh=180 [49] . 
The maximum energy of the shell correction is chosen to be Esh=5 MeV, and the widths 
of the shell are = az =5 MeV. These values are consistent with the magnitudes of 
the shell effects discussed for superheavy elements. In Fig. 9 we compare the calculated 
charge and mass distributions without and with the shell corrections (i.e., without and 
with the term AFaz)- We consider typical conditions in supernova matter where the 
production of superheavy nuclei is still possible, i.e., T =1 MeV, ps ~ O.OSpo, = 0.2. 
One can see that a pure liquid-drop description predict Gaussian-like mass and charge 
distributions of fragments which move to smaller values with decreasing density. When 
the shell constraint is imposed, the yields become essentially, by factor 2-5, larger in the 
vicinity of neutron [Ngh) and proton (Zgh) 'magic' numbers. Moreover, because of the 
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increased binding, these magic nuclei are abundantly produced in a rather broad region 
of density. For example, the neutron shell of Nsh=180 dominates clearly at densities of 
both 0.03po and O.OSpo- 

This result makes possible to discuss nucleosynthesis of superheavy elements at su- 
pernova conditions. As has been already pointed out in ref. [20], there is a chance that 
heavy nuclei could be produced at subnuclear densities, and then ejected into the space. 
At considered small temperature (1 MeV) the effect of secondary deexcitation will be 
minimal, only few nucleons will be lost. The fission channels for these neutron-rich nuclei 
may also be suppressed by the shell effects [50] and electron screening [51j. However, 
one can expect that these nuclei will fastly emit neutrons above the neutron drip-line, 
undergo abundant (3 decay and, possibly, a emission. We are planning to analyze all 
these processes in the forthcoming publications [52j. Here we mention only a possible 
scenario how such superheavy nuclei can be ejected into space. Since the synthesis of 
heavy and superheavy nuclei is only possible at rather high baryon density ps ~ O.OSpo it 
is most likely that such nuclei will not be ejected in the course of the supernova explosion. 
Instead, they will be accumulated at the surface of a newly produced neutron star. If 
this star is in a binary system with another neutron star, white dwarf, or even a black 
hole, there is a chance of their collision at later stages of evolution. Then a part of the 
stellar material will be ejected in space, while the other part may collapse into a black 
hole. One can also speculate about asymmetric explosions of supernovas, acceleration of 
nuclei by the neutrino wind, and starquakes which may provide this ejection (see, e.g., 
refs. [531 EH [55] )• The search for new mechanisms of nucleosynthesis is motivated by the 
fact that the traditional s- and r-processes have serious problems to explain synthesis of 
fissioning nuclei larger than lead [T3] . 

5. Thermodynamical characteristics of stellar matter. 

In this section we present general thermodynamical characteristics of stellar matter, 
such as energy density, pressure and entropy as functions of temperature T, baryon density 
Pb, electron fraction Y^, as well as the nuclear composition. In hydrodynamical simulations 
most important role is played by Equation of State (EOS), which connects pressure with 
the energy density. The dynamics of collapse and explosion depends essentially on the 
EOS [5j. 

5.1 Caloric curve, pressure and entropy 

One of the main inputs for dynamical simulations of supernova explosions is the total 
thermal energy deposited in the matter. In Fig. 10 we show so-called caloric curves, 
i.e., the thermal energy per nucleon of the matter as a function of temperature. For 
convenience, the actual energy per nucleon is shifted by the value of 16 MeV, which 
corresponds to the bulk binding energy of nuclear matter at ps = Po- Our calculations 
show that nuclear contributions dominate at high densities and low temperatures, where 
heavy nuclei survive. Due to the internal excitation of these nuclei according to the 
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compound nucleus law, E* ^ T"^, the caloric curve has a parabolic shape at ~ 0.1 — 
O.OOlpo- However, at low densities and high temperatures electrons and photons dominate, 
and the caloric curve behaves according to the Stefan-Boltzmann law E* ~ T^, i.e., 
the energy per nucleon grows very rapidly with temperature. One can also see that the 
nuclear contributions become nearly independent on baryon density at high temperatures. 
Namely, they approach the Boltzmann limit, E* ~ 1.5T, when nuclear matter disintegrate 
completely into nucleons. It is instructive to note that at high densities p ~ O.lpo and 
proton fractions Yp = ^ 0.4 which are typical for normal nuclei, the caloric curve is 
determined mainly by the nuclear species, and it reminds very much the caloric curves 
extracted from the multifragmentation reactions [2]. In this case, the excitation energy 
increases rapidly at temperatures T ^ 4 — 6 MeV, which correspond to the maximum 
in the heat capacity. This is a characteristic feature of the liquid-gas phase transition 

[Ml HZ]. 

Pressure is another important characteristic of the matter, which, in competition with 
the gravitational pressure, determines the structure and dynamics of the stellar system. 
In Figs. 11 and 12 we plot the pressure versus baryon density for different temperatures. 
There are several important features to be mentioned. First, the nuclear contributions are 
mainly important in the intermediate density region (10~^ — 10~^)po; where, depending 
on temperature, more and more free nucleons are present in the nuclear matter. In the 
case of full disintegration of nuclear species into nucleons, the nucleons may contribute 
up to 50% to the total pressure. Second, at higher densities the pressure is dominated 
by the relativistic electrons, since their chemical potential becomes very high. Third, 
at very low densities and high temperatures the radiation pressure dominate, which is 
proportional to T"^ and does not depend on baryon density. One can see that the nuclear 
pressure is higher in the case of a low electron fraction (Ye=0.2), because of a considerable 
abundance of free neutrons, even if large nuclei are present in the matter. In the modern 
hydrodynamical simulations of supernova explosions the shock stalls at densities around 
10~^ — 10~^po This happens partly because a significant fraction of shock energy is 
used for disintegration of infalling nuclei, from C to Fe. Therefore, survival of medium 
and heavy nuclei would contribute essentially to the revival of the shock. 

The entropy per baryon S/B, which is shown in Fig. 13 as a function of temper- 
ature, is an important characteristic of the exploding matter. One can notice that it 
correlates strongly with behavior of the caloric curve. At low temperatures and high 
densities the nuclear contribution to the entropy dominates. At high temperatures, the 
nuclei disintegrate into nucleons and the nuclear entropy depends only logarithmically on 
temperature and density according to the Boltzmann gas law. Usually, at entropy greater 
than 10 units per nucleon only nucleon gas without heavy nuclei is present in the system. 
The main contribution to the total entropy in this case is provided by the radiation and 
electron-positron pairs. This contribution does not depend on density and is proportional 
to T^. The total entropy has a jump across the shock, which can be explained in part by 
disintegration of heavy nuclei. Therefore, even small differences in nuclear properties in 
medium, in comparison with isolated nuclei, may lead to significant effects in the shock 
dynamics. 
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5.2 Adiabatic trajectories. 

Some important processes in stellar matter proceed with approximately constant en- 
tropy. For example, the collapse of a massive star before re-bounce is characterized by the 
entropy around one unit per nucleon. After propagation of the shock wave the entropy 
increases drastically. However, subsequent evolution of the matter is again close to isen- 
tropic, and this assumption is quite valid for nucleosynthesis. For this reason we consider 
adiabatic trajectories in the T — ps plane, and calculate fragment mass distributions and 
thermodynamical functions along these trajectories. 

In Fig. 14 we show adiabatic trajectories for several fixed S/B. This representation of 
the phase diagram onT — p plane is very convenient for understanding thermodynamical 
properties of stellar matter (compare also with Fig. 1). It is also possible to make a rough 
estimation of nuclear composition of the matter: As well known from previous studies 
at S/B=l many heavy nuclei exist in the matter, while at the entropy as high as 
S/B=20 the baryonic matter consist mainly of free nucleons and hadron resonances. Our 
calculations confirm this expectation. 

It is instructive to compare mass fractions of different nuclear species along the adi- 
abats for different thermodynamical conditions, in order to get an idea about evolution 
of nuclei during the whole collapse-explosion process. In Fig. 15 we show these fractions 
(similar to Fig. 6) for S/B=l, where the mass fraction of heavy nuclei is close to 1, and 
for S/B=8, where the nuclei undergoes deep disintegration. In the latter case the fraction 
of nuclei with A > 4 is essential only at high densities, and there are no large fragments 
at low densities. At S/B=8 the fraction of a particles has a very interesting behavior: 
It has a minimum at intermediate densities {ps ~ lO^^po); but then it increases with 
decreasing density, and completely dominates at low densities. For this fixed entropy a 
large binding energy of a particles is very important for the thermodynamical balance in 
the system. Therefore, a particles can be preferable instead of individual nucleons at the 
low density conditions. The production of a clusters may help to revive the shock wave 
by maintaining a sufficiently high temperature behind the shock. 

The mass distributions of nuclei evolve strongly along the isentropic trajectories. In 
Fig. 16 we demonstrate this evolution in stellar matter with a large electron fraction 
Yg = 0.4 for two cases corresponding to a low entropy S/B=l, and to a higher entropy 
S/B=4, where the contribution of heavy fragments is still essential. It is important that 
in the first case (S/B=l) at high baryon density (p^ = O.lpo) and temperature (T =3.39 
MeV), the distribution of heavy nuclei is centred around A ~130 and has a large width 
aA ~ 50. Therefore, heavy nuclei with mass number higher than 200, as well as very 
small clusters, coexist in the system. With decreasing density the distribution of nuclei 
shifts to smaller masses and becomes more narrow. Finally, the distribution moves into 
the iron region, where the binding energy is maximal. At higher entropy (S/B=4) we 
have an opposite situation. The masses of produced nuclei become larger with decreasing 
density. And at very low densities the nuclei reach finally the iron region. The reason for 
this interesting behavior is that the temperature drops essentially (from 10 to 1 MeV, as 
seen from Fig. 14) along the adiabatic trajectory. 
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All above discussed trends take also place in the case of isospin asymmetric matter 
with a small electron fraction, = 0.2, see Fig. 17. Besides of the expected effect of 
increasing number of free neutrons, we obtain here that more heavy nuclei are produced 
at the same densities. In the same time these nuclei are very neutron-rich. For example, at 
S/B=l and p = O.lpo) for Ye = 0.2 we have an average charge of big nuclei {Z)/A 0.28, 
in comparison with {Z)/A ^ 0.41 for Ye = 0.4. Because of formation of the unusual 
heavy nuclei, the temperature is lower for small electron fractions. One can conclude 
from analysis of Figs. 15, 16, and 17 that clusterization of nuclear matter at low densities 
should have important consequences for the explosion, at least via the energy balance. 
However, this effect will be even stronger if we include into consideration the modifications 
of the weak reaction rates caused by the clustering (see below). 

In Fig. 18 we demonstrate the total adiabatic pressure as function of baryon density. 
As expected, it increases with density, and with entropy. One can see a nearly linear 
relation between ln(P) and ln(p). For a fixed entropy this relation is usually expressed as 
the politropic equation 

p~p^''^ (19) 

where Tad is an effective adiabatic index. In Fig. 19 we show the behavior of Fad as a 
function of density for two values of entropy per baryon. At S/B=l, when heavy nuclei 
still exist, the adiabatic index is nearly constant and close to 4/3, which is expected for 
relativistic electron gas. At S/B=8 the Tgd coefficient shows more interesting behavior: 
it changes considerably and goes through the maximum. The maximal value of 1.5 is 
reached around 10~^ — 10~^po- This corresponds to a change of the pressure slope seen in 
Fig. 18. This effect is caused by a nearly complete disintegration of nuclei into nucleons, 
which takes place in this density region (see Fig. 15). The production of heavy nuclei at 
higher densities, and production of a particles at lower densities, lead to decreasing Tad 
at the both side of the maximum. Matter with Fad < 4/3 can not resist the gravity and, 
therefore, unstable with respect to gravitational collapse. On the other hand, matter with 
Tad > 4/3 can provide conditions for outward propagation of the shock wave during the 
supernova explosion. 

Finally, let us consider now the adiabatic sound velocity, = dP/dp\s, which plays an 
important role in hydrodynamical simulations. For example, when the collective velocity 
of matter exceeds a shock wave is generated. In Fig. 20 we present the sound velocity 
along different adiabates. In this case it can be obtained as = F^^ ■ P/ p. As expected, 
Cs increases with density. However, in the case of S/A=8, there is a peak around 10~^po) 
which is caused by the same physical reasons as the maximum of Tad, shown in Fig. 19. We 
note that a small sound velocity in comparison with the light velocity give a justification 
for using nonrelativistic formulas for baryonic and nuclear degree of freedom. 

6. Possible in-medium modification of nuclear properties. 

6.1 Reduction of symmetry energy. 
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As we have mentioned in section 3 multifragmentation reactions open a unique pos- 
sibility to investigate clusterization of nuclear matter at subnuclear densities. Recently, 
the symmetry energy of hot nuclei was extracted from experimental data [231 EH EHl EH], 
and it was demonstrated that the 7 coefficient, see Eq. ( IT2l) . is considerably reduced as 
compared with the values expected for cold isolated nuclei. This effect becomes stronger 
with increasing excitation energy. For example, in Fig. 21 we show the extracted val- 
ues of the 7 coefficient, which go significantly down with decreasing impact parameter 
b, i.e., for more and more central collisions. The empirical value of the 7 coefficient, 
approximately 25 MeV, was obtained for isolated nuclei from the liquid-drop description 
of their binding energies. As one can see, at high excitation energies it drops down to 
^ 15 MeV, and may even lower, as follows from the analysis of ref. [37j. As discussed 
in refs. [IHl EZl EHl EHl HOI IH] this change can be explained by the reduction of the 
fragment density, modification of the nuclear surface energy, and the influence of nuclear 
environment. 

Now we come back to the comparison of conditions which can be reached in multifrag- 
mentation reactions and in supernova explosions. Specifically, in Fig. 22 we demonstrate 
the similarity of fragment mass distributions for these two physical systems, calculated 
within the SMM and SMSM using the same description of hot fragments in dense en- 
vironment. The density and temperature of stellar matter are chosen close to the ones 
typical for multifragmentation reactions. The electron fraction {Y^ = 0.2) corresponds 
approximately to the deleptonization values obtained at this density. One can see that 
the evolution of mass distributions with excitation energy and temperature is qualita- 
tively similar for both cases (see also section 4.3 and refs. [21 EH]). The transition from 
the 'U-shaped' mass distribution to the exponential one, a characteristic feature of the 
liquid-gas phase transition, is well pronounced in both cases too. However, in the su- 
pernova environments much heavier and neutron-rich nuclei can be produced because of 
screening effect of surrounding electrons. In this figure we also demonstrate how impor- 
tant is to extract reliable information about the symmetry energy of hot nuclei. As one 
can see from mass yields at 3 MeV per nucleon in top panel, changing 7 coefficient from 
25 to 15 MeV has practically no effect on the mass distributions of fragments produced in 
nuclear reactions. As was shown in refs. [lOl EZ] , in nuclear multifragmentation reactions 
the most noticeable effect is that the isotope distributions become broader at smaller 7. 
However, the symmetry coefficient 7 has a dramatic infiuence on masses of nuclei pro- 
duced in supernova environments. One can see from Fig. 22 that much more heavy (and 
more neutron-rich) nuclei can be formed in this case. This effect makes very likely pro- 
duction of heavy and superheavy nuclei in supernova environments. In the following these 
hot nuclei should undergo de-excitation, and their decay products may either survive on 
the surface of a neutron star or be ejected into inter-stellar space. Also they can serve 
as seeds for subsequent r— process, as discussed in section 4.4. Therefore, studying the 
multifragmentation reactions in the laboratory is important for understanding how heavy 
elements were synthesized in the Universe. 

6.2 Electron capture rate. 
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Reduction of the symmetry energy of hot nuclei can be very important for weak pro- 
cesses. Here we consider a typical example related to deleptonization of matter, e.g., the 
electron capture by nuclei. One should bear in mind that the electron fraction is crucial 
for dynamics of stellar matter, since the electron pressure dominates at subnuclear den- 
sities. The calculation of the electron capture rate. Re, was carried out with the method 
suggested in ref. [58|. It is based on an independent particle model and assumes domi- 
nance of Gamow- Teller transitions. The electron chemical potential /ie and the reaction 
Q-value are the most important energy scales of the capture process. It is clear that the 
Q-value is directly related to the symmetry energy coefficient 7. A good approximation 
for the capture rate (per second) on an isolated nucleus is the expression [58j: 



where tg =6146s, Bg =4.6 represents a typical (Gamow- Teller plus forbidden) matrix 
element, ^ = (Q - 5E)/T, r] = {^e + Q - 5E)/T, 5E =2.5 MeV, and Fk are the rel- 
ativistic Fermi integrals of order k. It is instructive to normalize per nucleon this rate 
by taking into account the whole ensemble of heavy nuclei produced in stellar matter: 
(-Re) = H PAzRel Pb- Figure 23 demonstrates that the electron capture rate in stellar 
matter depends very essentially on the symmetry energy of nuclei. One can see that at 
relatively high densities pb ~ O.lpo the electron capture rate changes only by 20-50%, 
if we adopt the reduced symmetry energy coefficient 7 ^ 15 MeV. This is because a 
high electron chemical potential drives the reaction. However, at small densities (below 
10~^po); when heavy nuclei with large charge still exist (at least at low temperatures), 
the effect of reduced 7 is dramatic, of two-three orders of magnitude. We note, that 
at these relatively low densities and temperatures the nuclear chemical equilibrium may 
already be problematic [59] , although some authors keep using it in network calculations 
|60j . We believe that hot nuclei can interact with each other by neutron exchange in this 
case. This situation is similar to what we have at higher densities of nuclear matter in 
multifragmentation reactions. Therefore, the effect of reduction of the symmetry energy 
observed in multifragmentation may also take place in the supernova environments and 
be responsible for significant enhancement of weak reaction rates. 

Conclusions 

We have formulated a statistical approach (SMSM) designed to describe supernova 
matter at subnuclear densities. It may be applied for a broad variety of stellar processes, 
including the collapse of massive stars and supernova explosions, clusterization of nuclear 
matter in the crust of neutron stars, nuclear composition in merging binary stars, etc. 
The model includes the whole ensemble of nuclear species, as well as photons and leptons 
(e~, e"*", I/, v). The model fully accounts for the nuclear liquid-gas phase transition, which 
was previously under active investigation in nuclear reactions. In general, we emphasize 
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a close connection of the processes in stellar matter with multifragmentation reactions 
studied in laboratories. 

We have calculated main thermodynamical characteristics of the stellar matter under 
different assumption on the lepton fractions. Nuclear degrees of freedom contribute essen- 
tially to the energy and the entropy at high densities. Whereas, at low densities and high 
temperatures the photons and leptons contributions dominate. The nuclear contribution 
to pressure becomes essential only when nuclei completely dissociate into nucleons. Ac- 
cordingly, the adiabatic index increases considerably in this region. On the other hand, 
we have found that the a particle production at low densities and moderate entropies can 
be an important process, which should be correctly taken into account in the dynamical 
simulations of the explosion. 

The comparison with the Lattimer-Swesty model shows that thermodynamical quan- 
tities of the matter, e.g., pressure, are not very different in two models. The reason is 
that both models treat the leptons and photons in a similar way. Considerable differences 
appear in the yields of a particles and heavy nuclei. We believe that the SMSM pro- 
vides more realistic mass and charge distributions of hot nuclei, without any additional 
constraint on their sizes. 

As a result of our calculations, we especially emphasize the evolutionary nature of the 
mass and charge distributions of produced heavy fragments. These distributions carry im- 
portant information regarding the nuclear liquid-gas phase transition in the stellar matter. 
Also these nuclei participate in many processes which determine the energy deposition and 
dynamics of the collapse and explosion. Motivated by recent findings in nuclear multifrag- 
mentation reactions we have analyzed how possible in-medium modifications of nuclear 
properties, in particular, a reduction of the symmetry energy, can influence the fragment 
yields and weak processes. We have found that these effects can be very important, e.g., 
for the electron capture by nuclei, which is responsible for deleptonization of matter. At 
the same time, they can increase the yield of big neutron-rich nuclei. We have discussed 
new mechanisms of nucleosynthesis leading to the production of heavy and superheavy 
nuclei in supernova environments. In particular, the shell effects existing at relatively low 
temperatures may provide an additional enhancement factor for their formation. 
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PHASE DIAGRAM of NUCLEAR MATTER 
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Figure 1: Nuclear phase diagram on the 'temperature - baryon density' plane. Solid and dotted 
lines indicate boundaries of the liquid-gas coexistence region and the spinodal instability region. 
The shaded area corresponds to typical conditions for nuclear multifragmentation reactions. The 
dashed lines are isentropic trajectories characterized by constant entropy per baryon, S/B =1, 
2, and 4. 
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Figure 2: Estimated reaction rates for a typical nucleus, with mass number ^=60 and charge 
^=24, in stellar environment versus the density of free neutrons. Different lines correspond to 
different reaction types as indicated in the figure. Temperatures T (in MeV) are given at the 
lines. 




10 ® 10 ^ 10 10 , 10 
P/Pn 



' 10-' 10-^ 



Figure 3: Comparison of SMSM and the Lattimer-Swesty (LS) model [17J for stellar matter 
with electron fraction 1^=0.3 and temperatures T =1 and 2 MeV, as functions of the baryon 
density in units of the normal nuclear density po ~0.15 fm^3. The panels present the total 
entropy per nucleon, the pressure divided by density, the fractions of a-particles and heavy 
nuclei (^4 > 4). 
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Figure 4: Average fractions of electrons (top panels) and free neutrons Yn (bottom panels) 
versus baryon density in units of po- Left panels present results for the /3-equilibrated matter 
with T =1 and 3 MeV. Right panels correspond to the conserved lepton fractions 1^=0.1, 0.2 
and 0.3, and T =3 MeV. 
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Figure 5: Mass fractions of nuclear species for stellar matter at ¥^=0.2 as functions of tempera- 
ture. Solid lines are for heavy nuclei (^4 > 4), dashed lines - a-particles, dotted lines - neutrons, 
dot-dashed lines - protons. The results for baryon densities of 10~^, 10~^, 10~^, 10~^po are 
presented in the corresponding panels. 




Figure 6: The same as Fig. 5 but at Ye=OA. 
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Figure 7: Mass distributions, i.e., yields per nucleon, (bottom panel) and average charge to 
mass-number ratios (top panel) for nuclei produced at density 10~^po and ¥^=0.4, for temper- 
atures r=l, 2 and 3 MeV. 
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Figure 8: Fragment mass distributions (yields per nucleon) at T=l MeV and several densities 
in units of = lO'^po (see notations at the lines). Electron fraction is Ye=OA. 
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Figure 9: Mass (right panels) and charge (left panels) distributions of superheavy nuclei at 
subnuclear densities without (top panels) and with (bottom panels) shell corrections. The stellar 
matter has electron fraction Y(,=0.2 and temperature T=l MeV. Different lines correspond to 
different densities as indicated in the figure. 




Figure 10: Energy per nucleon as function of temperature in stellar environment with electron 
fractions 1^=0.2 (left) and 1^=0.4 (right), measured from the binding energy of normal nuclear 
matter (16 MeV per nucleon). Dashed lines show only contributions of nuclei. Solid lines give 
total energies including nuclear, electron and photon contributions. Baryon densities of 10~^, 
10~^, 10"*^, 10~^, 10~^, and lO^^po, correspond to the 6 lines from the top to the bottom in the 
both panels. 
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Figure 1 1 : Pressure as function of baryon density in units po (the phase diagram in the P — p 
plane) for Y(,=Q.2. Dashed hues are only nuelear contributions, solid lines are total pressures 
including also electron and photon contributions. The temperatures of 6, 4, 2, 1, and 0.6 MeV, 
are presented by different lines from the top to the bottom. 



10 



CO 

E 
> 

0) 



10 



-7 



10 



Y^=0.4 


1 1 ' 






total 
nuclear 


6 








y/>'' .. 4 

2 


T=6 MeV 




-"" 1..— . 


_i— — 








' 7=0.6 MeV 
















1 1 . 


1 1 


10-' 10"^ 


— 4 — 

10 10 
P/Po 


-2 -1 
10 10 1 



Figure 12: The same as Fig. 11 but for Ye=OA. 
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Figure 13: Entropy (per baryon) of stellar matter as function of temperature, for electron 
fractions Yg=0.2 (left) and Yg=0.4: (right). Dashed lines show only nuclear contributions, solid 
lines give total entropy including also electron and photon contributions. The densities of 10~^, 
10~^, 10~^, 10~^, 10~^, and 10~^pO) correspond to the 6 lines from the top to the bottom in the 
both panels. 




Figure 14: Isentropic trajectories on the 'temperature - baryon density' plane for different 
entropy per baryon values (S/B=l, 2, 4, 8, 20) indicated in the figure. Electron fractions are 
Ye=0.2 (top) and Ye=OA (bottom). 
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Figure 15: Mass fractions of nuclear species along isentropes with S/B= 1 (top) and 8 (bottom), 
for Ye=0.4:. Solid, dotted, dashed, and dot-dashed lines correspond to heavy nuclei {A > 4), 
neutrons, a-particles, and protons, respectively. 
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Figure 16: Mass distributions of nuclei along isentropic trajectories with S/B= 1 (top) and 4 
(bottom), for Ye=0.4:. Different lines correspond to specific baryon densities and temperatures 
(in MeV), as indicated in the figure. 
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Figure 17: The same as Fig. 16 but for Ye=0.2. 
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Figure 18: Adiabatic pressure as function of baryon density (in units of po), for Ye=0.2 (top) 
and Ye=0.4 (bottom). Different lines correspond to entropies per baryon S/B=l, 2, 4, 8, and 20 
units, as indicated in the figure. 
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Figure 19: Adiabatic index as function of baryon density at S/B= 1 (solid lines) and S/B= 
(dashed lines), for ¥^=0.2 (top) and Ye=OA (bottom). 
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Figure 20: Adiabatic sound velocity Cs (in units of light velocity) at S/B= 1 (solid lines) and 
S/B=8 (dashed lines), for Ye=0.2 (top) and Ye=OA (bottom). 
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Figure 21: The apparent symmetry energy coefficient 7 of hot nuclei, as extracted from 
multifragmentation of tin isotopes induced by ^^C beams with energy 300 and 600 MeV per 
nucleon, versus relative impact parameter b/bmax [37] . 
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Figure 22: Fragment mass distributions (yields per nucleon) in multifragmentation of Au 
sources (top panel) and in supernova environment at the electron fraction =0.2 and baryon 
density O.lpo (bottom panel). The calculations at excitation energies of 3, 5, and 8 MeV per 
nucleon (top), and different temperatures T (bottom), are shown by different curves. Effects of 
the reduced symmetry energy coefficient 7 are also demonstrated in both panels. 




Figure 23: Density dependence of electron-capture rates Re on hot nuclei in supernova environ- 
ment at different temperatures T and the electron fraction 1^=0.4. Solid and dashed lines show 
results for standard (25 MeV) and reduced (15 MeV) values of symmetry energy coefficients 7. 
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